Vlasov analysis of relaxation and meta-equilibrium 
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The Hamiltonian Mean-Field model (HMF), an inertial XY ferromagnet with infinite-range in- 
teractions, has been extensively studied in the last few years, especially due to its long-lived meta- 
equilibrium states, which exhibit a series of anomalies, such as, breakdown of ergodicity, anoma- 
lous diffusion, aging, and non-Maxwell velocity distributions. The most widely investigated meta- 
equilibrium states of the HMF arise from special (fully magnetized) initial conditions that evolve 
to a spatially homogeneous state with well defined macroscopic characteristics and whose lifetime 
increases with the system size, eventually reaching equilibrium. These meta-equilibrium states have 
been observed for specific energies close below the critical value 0.75, corresponding to a ferromag- 
netic phase transition, and disappear below a certain energy close to 0.68. In the thermodynamic 
limit, the ^i-space dynamics is governed by a Vlasov equation. For finite systems this is an approx- 
imation to the exact dynamics. However, it provides an explanation, for instance, for the violent 
initial relaxation and for the disappearance of the homogeneous states at energies below 0.68. 

I. INTRODUCTION 

Consider the one-dimensional Hamiltonian 

N N 

H = 2 + 2N £ t 1 - cos ^ - W ■ W 

i—l — \ 

It represents a lattice of classical spins with infinite-range interactions. Each spin rotates in a plane and is therefore 
described by an angle —tt < 8i < tt, and its conjugate angular momentum pi, with i = 1, . . . , N; the constant J is the 
interaction strength. Of course, one can also think of point particles of unitary mass moving on a circle. This model 
is known in the literature as mean-field XY-Hamiltonian (HMF) [H. 

The HMF has been extensively studied in the last few years (see Q for a review) . The reasons for such interest are 
various. From a general point of view, the HMF can be considered the simplest prototype for complex, long-range 
systems like galaxies and plasmas (in fact, the HMF is a descendant of the mass-sheet gravitational modeljj]). But the 
HMF is also interesting for its anomalies, be them model-specific or not. Especially worth of mention are the long-lived 
meta-equilibrium states (MESs) observed in the ferromagnetic HMF. These states exhibit breakdown of ergodicity, 
anomalous diffusion, and non-Maxwell velocity distributions, among other anomalies[3|, |J] (see also the contribution 
by A. Rapisarda et al. in this volume). It has been conjectured that it may be possible to give a thermodynamic 
description of these MESs by extending the standard statistical mechanics along the lines proposed by TsallisQ- 

The simplicity of the HMF makes possible a full analysis of its equilibrium statistical properties, either in the 
canonical[l| or microcanonical ensembles |6|. If interactions are attractive (J > 0), the system exhibits a ferromag- 
netic transition at the critical energy E c — 0.75JN. Here we will focus on the out-of-equilibrium behavior of the 
ferromagnetic HMF (J > 0), when the system is prepared in a fully magnetized configuration, at an energy close below 
E c , with uniformly distributed momenta ("water-bag" initial conditions). Under these initial conditions the system 
evolves to a spatially homogeneous state with well defined macroscopic characteristics and whose lifetime increases 
with the system size, eventually reaching equilibrium. Numerical experiments have shown the disappearance of the 
family of homogeneous MESs below a certain energy close to 0.68 JN. 
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II. EQUATIONS OF MOTION 

It is convenient to write the Hamiltonian JQl in the simplified form: 

i=l 

where we have introduced the magnetization per particle 

N 



y y^fj , with f 4 = (sin ^ , cos 6»i 



m N 

i=l 

and for simplicity we have taken J = 1. The equations of motion read 

• dH dH - 

^ = ^— = Pi, Pi = —57- = m-0;, 
dpi 96»i 

for i = 1, . . . ,N, with 0j = (cos^, — sinf?,). Without loss of generality, we can set the axes such that m x (t = 0) = 0. 
If, additionally, the distribution of momenta is symmetrical, then m x (t) — 0, Vi. In that case, the equations of 
motion become 

A dH . dH . A 

Notice that these equations can be seen as the equations for a pendulum with a time-dependent length. 



III. FIRST STAGE OF RELAXATION 



Fully magnetized states violently relax to a state of vanishing magnetization, within finite size corrections. The 
most elementary approach to describing the relaxation of m, from a given initial condition, is to perform a series 
expansion around t = 0, i.e., 

k>0 

In our case, the initial condition is such that m = 1 (with m, = 0), then one obtains the following coefficients for 
m y (t) 

c = 1 

C2 = -(p 2 )o 

c 4 = (/>o+4(p 2 ) 

ce = -((p 6 )o + 26(/)o + 16(p 2 )o + 18(p 2 )^ 



and c dd = 0, where averages are calculated with the initial distribution of momenta h(p). If h(p) at t = is 
symmetrical around p — 0, then m(t) is an even function of t. In particular, if the initial condition is water-bag, 
i.e., 9i = 0,Vi and additionally pi are uniformly distributed in the interval [—p ,p ], then (from Hamiltonian l|T|l. 
Po = Vfe, with e the energy per particle), one obtains 

m y (t) = l-et"+(^+ £ -y + ... . (2) 

The convergence of this series is very slow and, given that a general expression is not available, only the very short 
time of the relaxation can be described. 
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IV. VLASOV EQUATION 

On the other hand, the evolution equation of the reduced probability density function (PDF) in /i-space is formally 
equivalent to the Vlasov-Poisson systempj 

df df dvdf 

where V — — m • r(9) and m = J ddr(8) J dp f(9,p,t). If m = my, then 

at + p w ~ msm9 dp- = °' (4) 

with 

/7T POO 
d9 cosfl / dpf(9,p,t). (5) 
-7T J — OO 

The Vlasov equation (QJ can be cast in the form 

^ = [L + L x {t)]f, 

where Lq = —pde and L\(t) = m(t) sin9d p . We will consider states (for instance, with vanishing magnetization) for 
which the term Li(t) can be treated as a perturbation. It is convenient to switch to the interaction representation, 

i.e., to define /(*) = e^ Lot f{t), then 

dt yjJ ' 

where L\ = e" L °*Lie io *. 

The equation for the propagator U, such that f(t) = U(t)f(0), is dU /dt — L\XJ , therefore, 

U[t) = l+[ dt'Liit'pit'), 



and recursively, one has 

U{t) = 1+1 dt 1 L 1 (t 1 )+ I dhLx(ti) I dt 2 L 1 (t 2 ) + . . . . (6) 



The solution at order k of the Vlasov Eq. JIJ is 

f {k \9,p,t) = e-^m k \t)f(9,p,0), (7) 

where the index (k) indicates the order at which the expansion © is truncated. From here on, we will deal with 
continuous distributions, hence our treatment is valid in the thermodynamic limit. 



A. Lowest-order truncation 

At zeroth-order, the propagator is approximated by U ~ = 1. This is equivalent to neglecting the magnetiza- 
tion. Thus, if m = 0, the truncation is exact. 

For the initial distribution f(9,p,0) = g(9)h(p), where g(9) is uniform in [— tt,tt] (hence, m — 0) and h(p) is an 
arbitrary even function, both distributions remain unaltered in time, consistently with the numerical simulations in 
Fig. 7 of 0. In fact, if m = for any time, there are no forces to drive the system out of the macroscopic state. 

If the initial condition is f(9,p,0) = 8(9)h(p), where 8 is the Dirac delta function and h(p) an arbitrary even 
function of p (as in our case of interest), although m ^ 0, L\ is small (it is null at t = and remains small for later 
times), allowing a perturbative treatment. Then, we have 

oo 

f°\6,p,t) = e-P td °f(8,p,0) = h(p) 5(9 -pt) = —h(p) J2 z lk{e - pt) . 

k— — OG 
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Therefore, 

h<°>(p,t) = / def°\e,p,t) = h( P ), 



that is, at zeroth-order, the distribution of momenta, whatever it is, does not change in time. However the angular 
distribution does indeed change. For instance, in the particular case of the water-bag distribution, where h(p) is a 
uniform distribution in [—p ,Po], we obtain 



Ll 1 + 2j2cos(key- 

\ k>l 



Notice that the distribution of angles becomes uniform in the long time limit. It gets uniform through a mechanism 
of phase mixing, where particles do not interact (remember that magnetization has been neglected). From Eqs. (JSJ 
and JHJ), the zeroth-order magnetization is 

m<°>(t) = f d8 CO8 0gW(O,t) = (9) 

whose expansion in powers of time yields 

m(°)(t) = l-et 2 + ^ 2 i 4 + ... 
Observe that this expansion up to second-order coincides with the exact one, given by Eq. J2J), for any s. 

B. First-order truncation 

Now, recalling that = 1 + f Q dtiLi(ti), from J7J, at first-order, we have 

= fW(8,p 7 t)+e-P td » [ < ft 1 L 1 (t 1 )/(e,p,0), 



where L x {t) = cP td " \m(t) sin9 d p je-P tde . Then 

f (1) {0,p,t) = h(p)S(9-pt) + h'(p)S(9-pt) [ dhm(h) sin(6 + ptt - pt) 

Jo 

- h(p)S'(9-pt) / dh tim(ti) aia(e + pti - pt) . 
Jo 

After some algebra, for the case h(p) uniform in [— p ,p o ], we obtain 

m«(i) = m(°'(t) + f dt 1 1 ^(sm(p [t-t 1 }) + t -^sm( Po [t + t 1 })) (10) 
Jo 2 Po \ h + t J 

Substituting m(t) by m^°\t) one obtains the magnetization at first-order. Moreover, 

h (1) (p,t) = [ d6»/ (1) (6»,p,t) = h(p) + h'(p) [ dhmih) siniph) 
J -it Jo 

+ h(p) I dt\ txmiti) cos (pi i) . 
Jo 

We recall that, for the uniform distribution, h'(p) oc [8{p + p Q ) — 8{p — p )\- This explains why h(jp,t) presents two 
spikes at p = ±p„. 
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FIG. 1: Squared magnetization as a function of time. The initial state is fully magnetized with uniformly distributed momenta 
for e = 0.69. Symbols correspond to numerical simulations with TV = 1000, the distribution of momenta is regular. Dashed lines 
correspond to the zeroth-order approximation obtained from Eq. Full lines correspond to the first-order approximation 

given by Eq. 1101 



Fig. 1 shows the first stage of the relaxation of the magnetization. Numerical simulations were performed for 
N = 1000. Increasing the system size does not change the numerical curve in the time interval considered (t < 50). 
Of course, for longer times the curve becomes size dependent 0, Q The squared magnetization rapidly decreases from 
its initial value m 2 = 1 down to zero at t ~ 2. Then, it remains very close to zero up to t ~ 20. From then on, one 
observes bursts of small amplitude. Since Vlasov equation is exact in the thermodynamic limit, it describes the exact 
N = 1000 evolution up to time t ~ 50. The zero order approximation describes correctly m? vs t for a very short 
time (t ~ 0.4). The first-order approximation describes satisfactorily the violent initial relaxation (up to t ~ 2), but 
it does not reproduce the structure appearing later. Higher order corrections are required to describe that behavior. 
Extrapolation of numerical simulations 0, shows that m — > in the thermodynamic limit. This regime settles for 
times beyond the scope of our approximation. 



C. Equilibrium 

For completeness, let us discuss the distributions at thermal equilibrium [llj. If the system has already attained 
equilibrium, then dt = 0. Let also assume that the equilibrium distribution can be factorized, i.e., f(6,p) — g(9) h(p). 
Then, from @, 

Pff%) = msinW)^. (11) 

Assuming h(p) = Aexp(— /3p 2 /[2/i]), Eq. I|ll(l reduces to dg/d6 = — msmO g(6). Thus 

g (0) = Ce^ mcos0 , (12) 

with the normalization constant C — l/[2TrI (pm)], where lo is the modified Bessel function of zeroth-order. The 
equilibrium magnetization can be obtained from the consistency condition JSJ: 

m = d6 cos 9g{6) h{Pm) 



thus recovering the results of canonical calculations [jj. 
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D. Meta-equilibrium 

Although we have not found the long-time solution of Vlasov equation, starting from fully magnetized initial condi- 
tions, numerical simulations 3] indicate that in the thermodynamic limit the system tends to a spatially homogeneous 
state. We have seen in Sect. II V Al that. once reached a homogeneous state, the distribution of momenta, whatever it 
is, does not change in time. But, the question is whether the homogeneous solutions are stable or not under pertur- 
bations. On one hand, the Vlasov approach is a good approximation to the discrete dynamics, on the other finite-size 
effects may be the source of perturbations that may take the system out of a Vlasov steady state. Therefore, we 
will perform a stability test (valid in the thermodynamic limit) and discuss the results under the light of the discrete 
dynamics. 

There is the well known Landau analysis^ which concerns linear stability. A more powerful stability criterion for 
homogeneous equilibria has been proposed by Yamaguchi et al.|j|. This is a nonlinear criterion specific to the HMF. 
It states that f(p) is stable if and only if the quantity 




is positive (it is assumed that / is an even function of p). This condition is equivalent to the zero frequency case of 
Landau's recipe0,0- Yamaguchi et al.@ showed that a distribution which is spatially homogeneous and Gaussian 
in momentum becomes unstable below the transition energy e cr = 3/4 (see also 

in agreement with analytical 

and numerical results for finite N systems. They also showed that homogeneous states with zero-mean uniform f(p) 
are stable above e = 7/12 = 0.58... (see also [HQ)- I n the same spirit, it is instructive to analyze the stability of the 
family of g-Gaussian distributions 

f(p) ex e X p q (~ap 2 ) = [1 - a(l - q)p 2 } , (14) 

which allows to scan a wide spectrum of PDFs, from finite-support to power-law tailed ones, containing as particular 
cases the Gaussian (q = 1) and the water bag (q = — oo). In Eq. (|14fl . the normalization constant has been omitted 
and the parameter a > is related to the second moment (p 2 ), which is finite only for q < 5/3. In the homogeneous 
states of the HMF one has (p 2 ) = 2e — 1, as can be easily derived from Eq. (IJ. Then, the stability indicator I as a 
function of the energy for the g-exponential family reads 



3 — (7 , 
1 = 1- —, TT- r . (15) 

2(5-3g)(2e-l) V ' 



Therefore, stability occurs for energies above e q 



The stability diagram is exhibited in Fig. 2. It is easy to verify that one recovers the known stability thresholds 
for the uniform and Gaussian distributions. We remark that Eq. 1|16[) states that only finite-support distributions, 
corresponding to q < 1, are stable below e cr . This agrees with numerical studies in the meta-equilibrium regimes of 
the HMF. 

We have also shown recently 12] that a similar analysis can be performed for a very simple family of functions 
exhibiting the basic structure of the observed f(j>), basically, a uniform distribution plus cosine. Fitting of numerical 
distributions leads to points in parameter space that fall close to the boundary of Vlasov stability, and exit the 
stability region for energies below the limiting value e ~ 0.68. This result is confirmed when the stability criterion 
is applied to the discrete distributions arising from numerical simulations ^3 > although for the discrete dynamics the 
magnetization is not strictly zero. 

The stability index I is positive for energies above e ~ 0.68. The fact that the stability indicator becomes negative 
below e ~ 0.68 signals the disappearance of the homogeneous metastable phase at that energy. In fact, extrapolation 
of numerical simulations to the thermodynamic limit confirm this result. The present stability test only applies 
to homogeneous states. Strictly speaking, m = does not imply that the states are inhomogeneous. However, the 
sudden relaxation that leads to the present MESs mixes particles |4j in such a way that m — and spatial homogeneity 
are expected to be synonymous. Below e = 0.68, the measured distributions are evidently inhomogeneous (m ^ 0). 
In these cases, negative stability refers to hypothetical homogeneous states having the measured f(p). 
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q =-1.0 



FIG. 2: Stability diagram of the (/-Gaussian ansatz for the momentum PDFs. 



FINAL REMARKS 



We have seen that, although our approach is valid in the continuum limit, it gives useful hints on the finite size 
dynamics. Of course, it can not predict complex details of the discrete dynamics. However, the present approach gives 
information on the violent initial relaxation from fully magnetized states, for sufficiently large system. It also explains 
the disappearance of homogeneous MESs below a certain energy observed by extrapolation of numerical simulations 
to the thermodynamic limit. Moreover, the identification of MESs with Vlasov solutions is also consistent with the 
fact that when the thermodynamic limit is taken before the limit t — > oo, the system never relaxes to true equilibrium, 
remaining forever in a disordered state. 
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